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We consider the Josephson effect when two independent Bose superfluids are weakly connected, 
fn the presence of interparticle interaction and based on the calculations of the one-particle density 
matrix of the whole system, we find that the one-particle density matrix can be factorized which 
satisfies the general criterion of Bose superfiuid proposed by Penrose and Onsager. By introducing 
an effective order parameter for the whole system, our researches show that there is Josephson effect 
for two independent Bose superfluids. 

PACS: 67.40.-w; 67.40.Db; 32.80.Ys 



O ■ I. INTRODUCTION 

o 
u 

^ ■ In 1962, Brian Josephson predicted the Josephson effect^ by considering the quantum effect of two superconductors 
c/3 ' separated by a thin insulator. The Josephson effect is a natural quantum phenomenon for two coherently superposed 
macroscopic quantum objects. The subsequent experiments verified the Josephson effect and also lent support to the 
BCS theory about the physical mechanism of superconductors. The Josephson effect has been used in the invention 
of novel devices for extremely high-sensitivity measurements of currents, voltages, and magnetic fields. Up to date, 
there are still considerable interests in the fundamental physics and applications of this effect. Similarly to the case of 
superconductors, the ideal Josephson effect was also observed by investigating the fiow of superfiuid '^He through an 
Q ■ array of micro-apertures^. Recent experiments'^'^ also observed clear Josephson current and especially self-trapping 
' phenomenon^ for Bose condensates in dilute gases. 

About twenty years ago, P. W. Anderson raised a famous questional that when two initially separated superfiuids 
are connected, whether the two superfiuids would show a relative phase and therefore Josephson current. In the 
^ ordinary physical picture of Josephson effect, the initial quantum state is regarded as a coherent superposition of two 
macroscopic phase coherent quantum objects. In Anderson's question, however, before connecting two superfluids, 
, the two superfluids are completely independent. To answer Anderson's question, we shall consider the problem that 
. during the connection process or with the development of time after the connection process, whether the quantum 
state of the whole system will become a coherent superposition of two macroscopic quantum superfiuids. 

After the experimental realization of Bose- Einstein condensate in dilute gases^, there are a lot of interesting studies 
' about the splitting of a condensate and merging of two independent condensates, which relate closely to the above 
"p^. Anderson's question. The former question about the formation of a fragmented condensate during the splitting process 
d • of a condensate has been studied intensivelji^*2iiSiii. The quantum merging of two independent condensates (i.e., the 
^ I inverse process of the splitting) was also investigated recently by considering carefully the role of dissipation-^, and 
I ' by considering both the adiabatic and nonadiabatic mergin^S based on the well-known two mode approximation^^. 
' O Most recently, the splitting and the following merging of an elongated condensate^''' is considered by including the 
5 finite-temperature effect. On the experimental side, a continuous condensate source was created by periodically 
replenishing a condensate with new condensates^^. This striking experiment also gives strong motivation to study 
theoretically the merging of independent condensates. 

To describe the essential quantum feature of a Bose superfiuid, Penrose and Onsager— proposed the idea of off- 
diagonal long-range order (ODLRO) which gives the general criterion of a Bose superfiuid. The ODLRO plays a 
very important role in the description of the Bose superfiuid, especially because it has no classical analogi2*i&. For a 
5^ [ single superfiuid, if the one-particle density matrix (\E'''' (ri, i) ^ (r2, f)) can be factorized, i.e., {'^^ (ri, i) ^E" (r2, i)) = 
<&* (ri, t) <I> (r2, i), there is an ODLRO for the superfiuid, and the superfiuid can be regarded as a macroscopic quantum 
object which has stable spatial coherence propertjiii. For two initially independent superfiuids, here we investigate 
theoretically the dynamic process of the whole system when the barrier separating two superfiuids is decreased 
adiabatically so that the two initially separated and independent superfiuids are weakly connected. Based on the 
calculations of the one-particle density matrix of the whole system, it is found that there is an interaction-induced 
quantum merging process for initially independent superfiuids. After two initially independent superfiuids merge 
into a single condensate, it is shown that (ri, t) 5* (r2, t)) ~ (ri, i) $e (1*2 , t) with $e being an effective order 
parameter of the whole system. The effective order parameter has the property that it is a coherent superposition 
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FIG. 1: Shown in Fig. la is two independent Bose superfiuids in two separate tanks. After the two Bose superfluids are 
connected shown in Fig. lb, there is particle current between two tanks. 

of two macroscopic wave functions. For weakly interacting Bose superfluids, we find that the evolution of <i>e can be 
described very well by the ordinary Gross-Pitaevskii (GP) equation^^, and thus there is also Josephson current when 
two initially independent superfluids are weakly connected. 

The paper is organized as follows. In Section II we introduce the effective order parameter of the whole system 
when the one-particle density matrix is calculated for the general situation. In Section III we give the expression of 
the overall energy and the evolution equations based on the action principle. A brief summary and discussion is given 
in Section IV. 



II. EFFECTIVE ORDER PARAMETER FOR TWO INITIALLY INDEPENDENT BOSE SUPERFLUIDS 



For two independent Bose superfluids at zero temperature shown in Fig. la, the number of particles A'^i and in 
each of the two superfluids are fixed, and the corresponding quantum state is 

\m,N.)^^^_iairiair^\0), (1) 

where C„ is a normalization constant to assure {Ni, N2\Ni, N2) = 1. a[ (a^) is a creation operator which creates a 
particle described by the single-particle state 0i {(f)2) in the left (right) superfluid. In this initial quantum state, the 
quantum depletion is omitted. Thus, this initial quantum state is valid when a/l << 1 with a and I being respectively 
the scattering length and mean distance between particles. One should note that for two initially coherently-separated 
superfluids, the state is | A^) — (b^)^ |0) /^/W. with b'^ being a creation operator which creates a particle with the single- 
particle state [y'Ni(j}i + \/N2(f>2) /Vn with iV = A^i -I- A'2. As shown in Fig. lb, we consider the dynamic evolution 
process of the whole system when two initially independent superfluids are weakly connected. 

As shown in the following, after two independent Bose superfluids are connected, 4>i and 02 will overlap and 
especially become non-orthogonal in the presence of interparticle interaction. Thus, we consider the general case 
for J (j)i(f)2dV = ^ from the beginning. We first give the general expression of the one-particle density matrix 
(^''' (x, t) (y, t)) for this general case. Generally speaking, |C| << 1. However, after straightforward calculations, it 
is shown clearly that the nonzero value of ^ can play important role in the one-particle density matrix for A'^i |C| > 1 
and N2 Id > 1- 
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The operators ai and 0.2 can be written as ai = J ^(pldV and 02 = / ^'^j'^F' respectively. Here ^I* (x, t) is the field 
operator. By using the commutation relations of the field operators [S> (x, t) , 4* (y, t)] = and [i> (x, t) , ^I/^ (y, t)] = 
(5 (x — y), it is easy to get the commutation relation [oijOj] = C*- We see that Oi and are not commutative any 
more for J c/ii^j'^y being a nonzero value. 

It is well-known that the field operator should be expanded in terms of a complete and orthogonal basis set. 
Generally speaking, the field operator ^ (x, t) can be expanded as: 

$ (x, t) = ai0i (x, t) + k(t>'2 (x. t) + • • • , (2) 

where and 02 are orthogonal normalization wave functions. Assuming that (j)'2 = /? ('/'2 +Q:</'i), based on the 

r /2 

conditions / (pl(j)2dV = and / {(ji'^f dV = 1, we have = (l - ICI^) and a = - C*- Based on = / * ((/)^)* dV, 
we have 0,2 = k//3* + C^i. 

Because k and a\ are commutative, it is convenient to calculate the one-particle density matrix p (x, y, t) using the 
operators k and a\ . The exact expression of p (x, y, t) is 

p(x,y,i) = (iVi,7V2,t|$^(x,i)$(y,f)|iVi,iV2,i) 

= 71 -^t (x, (y, t) + 'y2e"^^ri (x, t) cj>2 (y, i) 

+72e-''^=</>^ (x, t) cj)^ (y, i) + 73.^; (x, i) </>2 (y, t) , (3) 

where the coefficients are 



71 = 



i=0 
72 = 51 



i=0 



73 



C2iV2!(iVi+*-l)!iVi(l-|C|') 


N2-Z „. 

|^|2. 


i\i\{Ni~l)\{N2~i)\ 




^ClN2\{N,+i)\(l-\C\'Y'~" 




i\{i + l)\{Nx-l)\{N2-i 


-1)! ' 


^ClN2\{N^+i)\{l-\Q\^Y'~" 





■e;' iUlNil {N2 - i - 1)\ ■ 



In addition, the normalization constant is determined by 



cl 



^^7V2!(^i+*)!(i-ici') ' icr 



i=0 



The phase factor ipc is determined by C = |CI 

The above onc-particlc density matrix is obtained based on the second quantization method. We have also proven 
that this one-particle density matrix is equal to the result calculated from the many-body wave function which satisfies 
the exchange symmetry of identical bosons. 

Introducing the ordinary order parameter $1 (x, = \/N\(j)\ (x, and $2 (x, = \/N24>2 (x, t), the one-particle 
density matrix can be naturally rewritten as 

p (X, y, t) = Pfac (x, y, t) + Prion (x, y, t) . (6) 

In the above equation, the factorable component p/ac(x, y, i) is 

p/„c(x,y,t) = $:(x,i)$e(y,i), (7) 
where we have introduced an effective order parameter 

$e (x, t) = (x, t) + ^/^e''^= $2 (x, t) . (8) 



The coefficients are 7^ = Ji/Ni, 72 = 72/\/^i^2 and 73 = 73/A^2- The non-factorable component Pnon (x, y,f) is 

Pnon (x, y,t)=v (e^^=<I>t (x, t) $2 (y, t) + e-^^^*; (x, t) $1 (y, t)) , (9) 
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FIG. 2: Proportion of the non-factorable component in the one-particle density matrix with different parameters. The relation 
between 71, 72, 73, and ^ for A^'i = N2 = 10'' is shown in Fig. 2a, while the relation between 71, 72, 73, and A''i = N2 for 
— is shown in Fig. 2b. For the case of Ni = N2, 71 = 73. For Ni(^ » 1, we have « 0, thus the one-particle density 
matrix can be factorized in this situation. 



where t] = ^73 — ^7(73 j . We see that the parameter rj sho-ws the proportion of the non- factorable component. If 

the coefficient r/ is approximate to 0, Pnon (x, y,t) can be omitted, and thus p(x, y, i) can be approximated as the 
factorable component pfac (x, y, t). 

For t"wo initially independent and ideal superfluids, because C = 0) — 0, based on the Schrodinger equation, it is 
easy to verify that ( (t) = at any further time. In this situation, p (x, y, t) — $* (x, t) $1 (y, t) + $2 ^2 (y, t)- 
Obviously, p (x, y, t) can not be factorized. Thus, the two superfluids can not be described by a single order parameter 
even there is an overlapping between two superfluids. 

When the interparticle interaction is taken into account, however, (t) can be a nonzero value. The nonzero value 
of C (t) physically arises from the interparticle interaction of the whole system. Although generally speaking, \C{t)\ 
is much smaller than 1 because (f'i4'2 is an oscillation function about the space coordinate, it is easy to show based 
on Eq. Q that a nonzero value of {t) plays a very important role in the density matrix for large Ni and -/V2. As a 
general consideration, shown in Fig. 2a is the relation between 7^, 72, 73, |?7| and C for Ni = N2 = , while shown in 
Fig. 2b is the relation between 7^, 72, 73, l^yl and A^i = N2 for (" — 10~^. We see that when Ni |C| ~ 1, the factorable 
component pfaci'^,y,t) gives significant contribution to p(x, y,t). In particular, when Ni |C| >> 1, one has \ri\ « 0, 
and thus the non-factorable component Pnon (x, y, t) can be omitted. In this situation, the one-particle density matrix 
can be factorized, and thus the whole system exhibits the property of ODLRO. According to the general criterion for 
Bose superfluid proposed by Penrose and Onsager, the factorability of the one-particle density matrix means that two 
initially independent superfluids have merged into a single superfluid which can be described by the effective order 
parameter <I>e (x, i). After two independent superfluids merge into a single superfluid, we see that the relative phase 
tpc emerges naturally during the interaction-induced coherence process. 



III. OVERALL ENERGY AND DYNAMIC EVOLUTION 
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We now turn to investigate the evolution equations of 0i and ^2- The overall energy of the whole system is 

E = y^dy(Ari,Ar2,i| (^-|^*tv2$ + Fext*^* + |$^$^$$^ |iVi,iV2,t), (10) 

where g = A-Kfr^a/m is the coupling constant with a being the interparticle scattering length. V^xt is the external 
potential of the system. 

After straightforward calculations, the overall energy of the whole system is 

E = Ekin + Epot + Eint, (11) 

where the kinetic energy Ekin is given by 



■ ^"^^ + ^^"^^ • V02) , (12) 

while the potential energy is given by 

= J Vextp{r,r,t)dV. (13) 
In addition, the interaction energy Eint of the whole system is given by 
Eint = J {N,,N2,t\^%^^^%%\Ni,N2,t)dV 

= I y" dV [hi\<t>if + h2\<t>2f + hs\<Plf \hf 

+Re (hi \<Pif cj>\cl>2e"^' + {<Plf ^e''"^'' + he |02|' '^^<^2e'^»)] , (14) 

where the coefficients are 

, ^ ciN2i{N, + i-2y.m{N,-i) ,.2\^=-,^|2i 



i=0 

N2-2 



i\i\{Ni-2)\{N2-i)\ 
ClN2\{Ni+i)\ / ,^,2\^2 



- 1^ z!z!7Vi!(iV2-.-2)!i^"l^l ) '^1' 

AClN2\{N,+i-l)\N, |.2\^=-*-l ^|2i 

= 1. ^^^^i^N,-X)\{N2-^-\)\V-\^\) '^1' 

. 4fi^,jV2!(iVl+^-l)!iVl ,^,2X^-^-1 ^,2i+i 

.!(. + l)!(iVi-2)!(iV2-z-l)!l^~'^'i ^ 



i=0 



h - Sr 2ClN2\{N,+i)\ / _ |^|2\^--2 

''^ - 2^,\(i + 2)\(N,-2)\{N2-i-2)\V 1^1; 



z!(i + 2)!(7Vi-2)!(7V2-i-2) 



i=0 

^'"^ 4C^Af2!(iVi+i)! / „xJV2-i-2 
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It is well known that the action principle is quite useful to derive the time-dependent GP equation for a single Bose 
superfluid. By using the ordinary action principle and the above overall energy, one can get the following coupled 
evolution equations for and (j)2- 

d(j)i ^ J_SE_ 
dt Ni Scbl ' 

^H'4l . (16) 

where 5E/5(t)\ and 6E/6(j)2 are functional derivatives. 

Although the coupled evolution equations lfTB|) are quite complex, for the case of A^i |C| >> 1, |CI >> 1 and 
Ni ~ 7V2, there is a very concise approximate evolution equation. When these conditions are satisfied, the overall 
energy of the whole system can be approximated very well as 

where 

ELn = 1^ / ^'^e • V<f ed^, (18) 



^pot — 



I 



(19) 



E' 



(20) 



For the cases of Ni |C| >> 1, N2 |C| >> 1 and Ni ^ N2, it is easy to verify that i5^j„ « Ekm and E' 



pot 



Epot based 



on the analogous analyses about the effective order parameter that p (x, y, t) w $* (x, t) <i>e (y, t) in this situation. For 



iVilCI >> 1, iV2|CI » 1 andTVi 
expanded as: 



N2, one can also prove the result of £'j'„( w -Emt. Based on Eq. (|20|l . E'^^^. can be 



TP' 



+Re (/34 



t>2e 



/32i</'2r+/?3i</'iri02i 

/35(0I)''^2e''^= 



lip 



■P6\(l>2\ <^t<^2e 



life 



(21) 



where 



f3i 


= ll 




= ll 




= 47173, 


Pa 


= 471^7173, 


P5 


= 27173, 




= 473^7173. 



(22) 

To compare with the exact expression of the interaction energy given by Eq. (|14|l . Fig. 3 shows the relation between 
h,/f3, (i = 1, • • • ,6) and C for Ni = N2 = 10^. It is shown clearly that for A^i |C| >> 1 and N2 |C| >> 1, hJP, « 1, 
and thus w Ei„t- 

Therefore, the overall energy can be approximated well as 



E^-^ J KWH,dV + J Ve^t\<fe\^dV + ^ J dV\<^e\\ (23) 

where $e (x, t) is the effective order parameter given by Eq. |(HJ). Based on this approximate energy expression, it is 
easy to get the following evolution equation for <I>e (x, t) : 

Qg) J;2 

in-g^^- — VHe + Ve.t^e+9\<i>e\''<i>e- (24) 
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FIG. 3: Shown is the relation between hi/ Pi (i = l,---,6) and C for A^i = iV2 = 10^. We see that for iVi \C\ » 1 and 
iVa Id » 1, h,/P^ ^ 1, which 



It is quite interesting to note that the above equation has the same form as the standard GP equationiS. 

Because there is no exact analytic solution for Eq. H15|l . here we give the solution based on numerical calculations. 
We consider the one-dimensional dynamic process when two initially independent condensates in dilute atomic gases 
are weakly connected through adiabatically decreasing the height of the barrier separating two condensates. The 
time-dependent double-well potential is assumed as 

(x, t) = \^muJix^ -h U (t) e-^'/^' , (25) 

where U (t) — [/qg^^* -f- Ui. The first term is the external potential due to a magnetic trap or an optical trap, while 
the second term is the central barrier due to a laser beam separating two condensates. When the central barrier 
separating two condensates is sufficiently high so that there is no tunneling current, one may prepare two completely 
independent (rather than coherently separated) condensates by directly cooling the dilute gases in the double-well 
trap. In the last ten years, the remarkable experimental advances^ on Bose condensate in dilute gases make it be 
quite promising to experimentally test the theoretical predication in this work. 

In the numerical calculations, we introduce the dimensionless parameters xq — x/lx^ t ~ Eit/h, g' ~ g/Eil^. Here 
Ix = ^/hjmujx and Ei = f? jlmll.. The particle number is A^i = 2 x 10'* and = 10^. In addition, t/p = 400i?;, 
\J\ — 2QEi, 9 = 2Ei/h, Ax = and g' Ni = 50. For these parameters, the tunneling effect can be omitted for two 
independent condensates at the initial time. When the central barrier due to the laser beam is decreased, there is 
particle current between two condensates. 

In the numerical calculations, first we get the ground state wave functions 0i and 02 at the initial time. Then, by 
numerically calculating the coupled equations we obtain the evolution of ^ {t) based on the numerical results of 
and 02- From C,{t), we give the evolution of \rj\ in Fig. 4. We have also confirmed in the numerical calculations 
that, for g' — 0, the numerical result of \C,\ can be regarded as zero because it is smaller than 10^^". This shows 
clearly that the nonzero value of |C| physically arises from the interatomic interaction, rather than the error in the 
numerical calculations. In fact, if we assume that C {t) is always zero with the development of time, it is easy to show 
that this is an inconsistent assumption because with this assumption we prove both in the numerical calculations and 
analytic analysis that |C| will increase from zero after the overlapping between two condensates. In the inset of Fig. 4, 
we give the numerical result of A^i {t) /N which shows clearly the Josephson effect. Our numerical results show that 
{Ni (t) + N2 {t)) /N is always equal to 1 with an error below 10~^^ which confirms further our numerical calculations. 

Although the Josephson effect shown in Fig. 4 is the numerical results of Eq. 1)16(1 . one can understand easily this 
effect with the effective order parameter. Due to the quantum merging process, when the central barrier separating 
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FIG. 4: Quantum merging process shown through the evolution of I77I and the Josephson effect. The evolution of Ni/N shows 
clearly the Josephson effect. Due to the particle current between two initially independent condensates, as shown in the inset, 
|»7| decreases significantly with the development of time. At t > 3, {rjl « 0, thus two initially independent condensates have 
merged into a single condensate. 



two independent condensates is decreased, two initially independent condensates will merge into a single condensate 
described by the effective order parameter <l>e whose evolution is determined by GP equation (|24|l . Thus, it is natural 
that there is coherent Josephson current for two initially independent interacting condensates. 



IV. SUMMARY AND DISCUSSION 



In summary, we investigate the dynamic process of the whole system when the barrier separating two initially 
independent Bose superfluids is adiabatically decreased so that the two initially separated and independent superfluids 
are weakly connected. When the interparticle interaction is considered, we show that there is an effective order 
parameter for the whole system under appropriate condition. Compared with previous interesting studiesiSii^ii^ about 
quantum merging, in the present work, we consider this problem by stressing the non-orthogonal property between 
the wave functions for different condensates after their overlapping and in the presence of interparticle interaction. 
In particular, it is found that the effective order parameter satisfies the ordinary Gross-Pitacvskii equation, which 
means that there is also Josephson effect for two initially independent Bose superfluids. This result for the effective 
order parameter makes our theory can be tested and applied directly in future experiments about quantum merging 
process, such as the experiment about Josephson effect for independent Bose superfluids. 

We stress here again that, in our theory, the quantum depletion originating from the elementary excitations at 
zero temperature is omitted in the initial state QJ. Based on the Bogoliubov theory of the elementary excitations^^, 

the number of particles due to the elementary excitations is of the order of {a/Tj and thus the quantum depletion 
is negligible for Bose condensate in dilute gases. There is another reason why the elementary excitations can be 
omitted in the effective order parameter. A simple analysis shows that and {(l>2\4>k) (here (pk is the normalized 

wave function of the elementary excitations) are of the order of |^| e"^'^"^'^/^-' with k and L being respectively the 
wave number of the elementary excitations and spatial size of the Bose superfluid. This exponential decay of ((/>i|0fc) 
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and {4>2\<t'k) originates from the integral where there is spatially oscillating phase factor in the wave functions of 
the elementary excitations and Bose superfluid. Thus, even there are elementary excitations, its contribution to the 
effective order parameter is negligible. Although the present theory can not give quantitative predication for liquid 
superfluid of *He because the quantum depletion is very important, we believe that the quantum merging process 
means that after two separate tanks are connected by a pipe, two supcrfluids of "'He can merge into a single superfluid. 

For the quantum state given by Eq. (1), with the development of time, the wave functions and (/)2 are no more 
orthogonal in the presence of interparticle interaction. If we use the orthogonal bases 01 and (p^ , the quantum state of 
the whole system is |A^i,A^2) ~ ^\)^'^{k/ [3* -\-Cjii)^'^ |0). We see that the number of particles in the orthogonal modes 
(j)i and 02 are no more definite. This quantum state becomes a superposition of different number of particles in the 
orthogonal modes (pi and (j)^ ■ This implies strongly that in some sense our theory is equivalent to the Gutzwiller type 
approach^ where the coupling between different modes leads to the coherent transfer of particles between different 
orthogonal modes. 

It is straightforward to generalize the present work to the quantum merging process of several independent Bose 
superfluids. It is possible that this quantum merging process contributes to our understanding of the Bose condensation 
process. During the evaporative cooling process, firstly there would be a series of independent subcondensates formed 
from the thermal cloud. Due to the quantum merging process, these independent subcondensates will overlap and 
finally merge into a single condensate with well spatial coherence property. During the evaporative cooling process, 
due to the thermal equilibrium of the whole system, the thermal atoms continuously jump into the ground state. The 
quantum merging process has the effect that the atoms in the ground sate merge into a single condensate which has 
well spatial and phase coherence property. 
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